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Summary. Effective models are derived from the ferromagnetic Kondo lattice 
model with classical corespins, which greatly reduce the numerical effort. Results for 
these models are presented. They indicate that double exchange gives the correct 
order of magnitude and the correct doping dependence of the Curie temperature. 
Furthermore, we find that the jump in the particle density previously interpreted as 
phase separation is rather explained by ferromagnetic polarons. 

Manganites [1] are often described by the ferromagnetic Kondo lattice 
model, which is considered to explain some of their features, e.g., the transition 
from antiferromagnetic to ferromagnetic order with doping [2] . 

The application of the model is motivated by the fact, that crystal field 
splitting divides the five d-orbitals into two e g and three t^g orbitals, where 
the latter are energetically favored in the case of manganites. All three ti g 
orbitals are singly occupied and rather localized. Due to a strong Hund's 
rule coupling, these electrons are aligned in parallel and form a core spin 
with length S = 3/2. The filling of the e g orbitals is determined by doping 
and these electrons can hop from one Mn ion to the next via the intermediate 
oxygen. Hund's rule coupling leads to a ferromagnetic interaction between the 
itinerant e g electrons and the ti g core spin. The core spins interact through 
super exchange leading to a weak antiferromagnetic coupling between them. 

In this chapter, we derive effective models for the ferromagnetic Kondo 
lattice model and introduce suitable Markov chain Monte Carlo (MC) algo- 
rithms. The presented results, were not obtainable by simple analytic con- 
siderations, are partly found by this MC method and partly by use of the 
Wang-Landau algorithm [3] . 
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1 Model Hamiltonian 

The spin and charge degrees of freedom in manganites can be described by the 
ferromagnetic Kondo lattice model with two orbitals (x 2 — y 2 and 3z 2 — r 2 ): 

H = - ^2 tia -3l3 4aa C jf)a ~ ^ <T ia ■ S t + J' ^ S t - Sj , (1) 

i,j,ct,/3,cr iot <ij> 

where c\ aa (c iaa ) creates (annihilates) an electron with spin a in orbital a at 
site i, <Ti a denotes electron spin in orbital a and Si the core spin. The first 
term of the Hamiltonian describes the hopping between the nearest neighbor 
sites; the hopping strength Uajp depends on the involved orbitals and the 
direction. As matrices in the orbital indices a, f3 — 1(2), corresponding to the 
x 2 — y 2 (3z 2 — r 2 ) orbitals (see e.g. [2]), this hopping reads 

t t (° \ t -t ( 3/4 TV3/4 ^I (2) 

- t y Ql j, ti,i+«/D-t ^ Tx /3 /4 x/4 ) ■ w 

The overall hopping strength is t, which will be used as unit of energy, by 
setting t — 1. 

The second term contains the ferromagnetic interaction between the elec- 
trons and the core spins and the third term is the AFM superexchange of the 
core spins. 

The S = 3/2 core spin can approximately be treated as a classical spin, 
which corresponds to the limit S* ^ oo [4, 5]. It is then replaced by a vector 
of unit length and the factor 3/2 is incorporated into Jjy. This approximation 
simplifies calculations enormously and should not lead to much difference from 
the quantum case except possibly for very low temperatures T « [5-9] . 

For large Jh, the electronic density of states of the Hamiltonian (1) is split 
into the lower and upper Kondo band, separated by approximately Jh- The 
e g electrons move mostly parallel to the core spins in the lower band, and 
anti-parallel in the upper band. 

In order to derive effective low-energy models for Jh >• t, J', we change 
the quantization axis for the electron spin from the global quantization axis 
(e.g. the z-direction) to a local quantization axis, namely the direction of the 
local t2 g core spin. Spin up (down) then means that the e g electron spin is 
parallel (antiparallel) to the core spin. The Hund's rule term Jh o- ia ■ Si 
becomes Jh (jii a [ — ^iaj), with the factor 1/2 coming from the electron spin 
also absorbed into Jh- While the e 9 -spin is preserved in global quantization, 
this is no longer the case in the local quantization. An up electron at site 
i can therefore become a down electron at site j, which is denoted by the 
superscript for the hopping strength. Furthermore, the hopping now depends 
on the core spins: 

Cj/S = W/3%T ( 3 ) 

The first factor is the orbital-dependent hopping strength (2) and the second 
factor contains the relative orientation of the core spins: 
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ul'J = CiCj + SiSje™^-^ = cos(i?y/2) e**« (4) 
u*'-" = a { Cl s 3 e- 1 ^! - CjSi e' 117 ^) = sm(^/2) ^ (5) 

with the abbreviations Cj — coa(9j/2) and sj = sm(0j/2) and the restriction 
< 9j < 7r, where #,,0, are the polar coordinates for core spin Si. These 
factors depend on the relative angle z?^ of the core spins Si and Sj and 
on some complex phases ipij and Xij- With a shift of the chemical potential 
/x — > /x — Jh, the Hamiltonian (1) in local spin-quantization reads: 

6 = - E C^4 CT s^+ 2J «E^i+ J 'E^-^- ( fi ) 

i : j,a : l3,fT,(7' ia <ij> 

This is still the same Hamiltonian as (1) without any approximation besides 
the use of classical core spins. 

1.1 Effective Spinless Fermions 

Most of the experimental results on manganites and all of the theoretical work 
presented here concerns electron densities < n e \ < 1, i. e. predominantly 
the lower Kondo band. As Jh is much larger than the hopping t and the AFM 
superexchange J', one can simplify the model by a separation of energy scales 
[10]. As a first approximation, one can take Jh — > oo and thereby leave out 
the configurations with e g electrons antiparallel to the core spins completely. 
This approximation is widely used [11,12], but misses some important effects 
discussed in Sect. 3.1. However, if one treats these configurations in second 
order perturbation theory [13,14], almost perfect agreement to the original 
Kondo lattice model is obtained without any additional numerical effort [15]. 
This approach is similar to the derivation of the t—J model from the Hubbard 
model, while the Jh — > oo method corresponds to U — > oo for the Hubbard 
model. 

In this effective model, the dynamical degrees of freedom are the low energy 
states with the e s -spins parallel (i. e. up in local quantization) to the t 2 g-spins. 
The virtual excitations meditated by the hopping matrix are configurations 
where one e g electron is antiparallel (down): 

(ia |) — > (j/3 I) — > (i'a 1 |). As the low energy states contain only up electrons, 
this lead to an effective spinless fermion Hamiltonian: 

i,j,a,/3 i,a,a' / /»\ 

_ E I E 2J H J ) C i>a> C ia + J' E &j ■ 

[ijti'],a,a\j,0 ' / <ij> 

The first term of this Hamiltonian contains the kinetic energy of the electrons 
moving in the lower Kondo band. As iH is largest for parallel core spins, 
this term favors ferromagnetism. The second term describes electrons that 
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get excited into the upper Kondo band and then hop back to the original 
site. It yields a density dependent antifcrromagnctic interaction between the 
core spins. The third term is a 'thrcc-sitc-term' of minor influence and will be 
neglected [15]. On the other hand, its inclusion does not increase the numerical 
effort. 

The reduction of the Hilbert space achieved by this effective model is the 
same as for the Jh — ► oo limit, and finite Jh can thus be treated with the 
same numerical effort. 

1.2 Uniform Hopping Approach 

A significant further simplification is the uniform hopping approximation pro- 
posed by van den Brink and Khomskii [16]. This approximation replaces the 
different angles of neighboring core spins by a mean value. In order to treat 
anisotropics, two different angles are chosen, 9 Z — Si - Si± z in z-direction and 
9 xy = Si ■ Si± x = Si ■ Si± y within the xy-plane. These should not be confused 
with the polar angle of an individual core spin It is assumed that the rela- 
tive orientation is the same between all nearest neighbor pairs. The hopping 
matrix therefore becomes translationally invariant. Spin configurations that 
are still treated exactly include, among others, ferro- and antiferromagnetism 
and spin canted states. 

The impact of the core spins on the hopping simplifies to 

<' ct = cos(|h Uz , <>- CT = sin(|) = v ^^2 (8) 

in z-direction and analogously in x/y-direction. Likewise, the inner product 
of the t2g spins entering the superexchange term can be expressed by 

Si-Si+ a = cos^ = 2«2_i. (9) 

The energy of this model can easily be evaluated, especially in the ther- 
modynamic limit and the ground state can be obtained by minimizing the 
energy with respect to 9 Z and 6 xy . 

For a one orbital model in one dimension with periodic boundary condi- 
tions, the Hamiltonian simplifies to 

H = -u^ 4c, ]T cjq + J'L(2ul l) , (10) 

(ij) * 

This Hamiltonian yields a shifted tight-binding band structure 

e k = -2ucos(/c) - (l-u 2 z )/J H (11) 

with a band width of 4 m 2 , which has its maximum for ferromagnetic core 
spins and vanishes for antiferromagnetic order. 



Aspects of the FM Kondo Model 



5 



Similar calculations can be done for three dimensions with both orbitals 
and a ground state phase diagram can thus be obtained, see [15, 16]. 

This UHA approach was extended to finite temperatures in [17]. To intro- 
duce this method, we proceed as follows: For a given core spin configuration 
S, characterized by the set of angles {9i, we define the average u-value 

^) = ]v-E u l/( 5 )- ( 12 ) 

P (ij) 

Here N p is the number of n.n. pairs (ij). In the ESF Hamiltonian (7), u\j 
is then replaced by u(S). Besides u\j the Hamiltonian depends on \u^ a \ 2 
and on Si ■ Sj, which correspond to sin (i?{j/2) and cost?^, respectively. As a 
further approximation, these terms are replaced by l~u 2 (S) and 2u 2 (S) — 1 
respectively, which leads to the one-orbital UHA Hamiltonian 

H = -u ]T cj Cj - ^ ]T z lC \ Cl + J'N p (2u 2 1) , (13) 

(ij) 1 

with Zi being the number of nearest neighbors for site i. This Hamiltonian 
defines the Boltzmann factor for the spin configuration S. In order to calculate 
thcrmodynamical expectation values, one still needs to calculate the density 
of states -T(u), i.e. the number of spin configurations with the same average 
value for u. It can be calculated exactly in one dimension and by use of the 
Wang-Landau algorithm [3] in higher dimensions. 

Once the density of states r(u) has been obtained, observables can be 
obtained for any temperature, much larger lattices can be treated and a 3D 
phase diagram for finite temperatures can be obtained. The numerical effort 
is reduced from an integration over the L-dimcnsional space of the core spin 
configurations to an integral over the one-dimensional unit interval for u. 



2 Monte Carlo Algorithm 

The algorithm used to simulate the Kondo lattice model and the effective 
spinless fermion model in the grand canonical ensemble is the one proposed 
in [18]. For each core spin configuration, the resulting Hamiltonian for the e g 
electrons is a one-particle problem. The statistical weight for the core-spin 
configuration S in the grand canonical ensemble is the starting point of grand 
canonical Markov chain Monte Carlo simulations: 

t -P(H(S)-ij,N) 

«'<^> - ' (14) 

It is calculated from the eigenvalues of H(S) by use of free fermion formulae, 
which is denoted by the trace over the fermionic degrees of freedom tr c . 
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As some particle numbers are not stable in the grand canonical ensemble, 
we developed a canonical algorithm. An exact approach would mean calcu- 
lating the Boltzmann weight for every possible distribution of electrons 
on L eigenvalues and summing over these contributions. This is numerically 
too demanding. But for low temperatures, only very few of these distributions 
actually contribute to the partition function. They can be obtained by filling 
No < iVgi electrons into the A^o lowest eigenenergies and distributing only the 
remaining N e \ — Nq s=a 5 electrons on the states around the Fermi-energy. The 
canonical weight then becomes: 

W (S\N el) = , (15) 

where V denotes these restricted permutations. In order to decrease autocor- 
relations, particle fluctuations within a set of 3 to 5 densities were allowed. 
The amount of core-spin rotations was small for most updates in order to en- 
sure high acceptance, but occasionally a complete spin flip was proposed. In 
one dimension, whole sections of the chain were rotated at once. 50 to several 
hundreds of sweeps were skipped between measurements. This ensured statis- 
tical independence for the ID calculations, in 2D, remaining autocorrelations 
were treated by autocorrelation analysis. 

While observables depending only on Z are independent of the spin quanti- 
zation (global/local), care must be taken when evaluating e.g. the one particle 
Greens function, which in global quantization can be written as 

« 4 »«= / V[S] w(S\n)u)] (S) « c t ; cj »£ , (16) 

where <C c i ; is the Green's function in local spin quantization. It can be 

expressed in terms of the one-particle eigenvalues and the corresponding 
eigenvectors ip^ of the Hamiltonian H(S): 

<< 4 ' -' >>w ^ u-{eW-n)+iO+ 1 j 

It should be pointed out that the one-particle density of states (DOS) is 
identical in global and local quantization; for details see [17]. 



3 Results 

In the first subsection, we will demonstrate the validity of the simplified ef- 
fective models. We will then use the uniform hopping approach (UHA) in 
Sect. 3.2 to determine the Curie temperature of the one-orbital model in 
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three dimensions. In Sect. 3.3, we will present results for the one- and two- 
dimensional model and we will show that they do not indicate phase separa- 
tion, but rather ferromagnetic polarons. A phase diagram for the 2D model 
is given in Sect. 3.4. 



3.1 Validity of the ESF Model and the UHA in one Dimension 

In this section, we will present results from unbiased Monte Carlo Simulations 
for the ESF model in one dimension with one orbital and open boundary 
conditions and we will compare them to results for the full ferromagnetic 
Kondo lattice model and the Jh — ► oo approximation. Simulations were done 
for L = 20 sites, = 50, J' = 0.02 and J H varying from 4 to 10. 




Fig. 1. Spin structure factor for the t-2g spins at n = 1 (inset: n ~ 0.75) for 
f3 = 50, J' = 0.02, L = 20 and different values of Jh- Circles: spinless fermion model 
(7); crosses: DE model (6). In the limit Jh —* oo (dashed line) the intensity of the 
AFM peak is considerably smaller than for finite Jh- From Ref. [15]. 



Figure 1 shows the core-spin structure factor for the three models. For 
electron density n w 0.75 (inset), the ti g correlations are ferromagnetic, driven 
by the kinetic energy of the e g electrons. For n = 1, the lower Kondo band 
is completely filled, no hopping is possible and the kinetic energy therefore 
vanishes. Excitations into the upper Kondo band (virtual for the ESF), which 
are favored by antifcrromagnctism, then dominate the energy. 

For both densities and even for moderate Jh = 4, the ESF model (7) 
and the original Kondo Model (6) produce virtually identical results. The 
Jh — * oo model on the other hand does not reproduce the antifcrromagnctic 
correlations for the completely filled lower Kondo band correctly, because 
the virtual excitations are missing from this model, and it also overestimates 
the ferromagnetic correlations at n w 0.75. The AFM effect coming from the 
virtual excitations can be described by a density dependent effective parameter 
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J e g- — J' + 1/(2 J H ) at n = 1 and it is generally much stronger than the small 
superexchange J' also favoring AFM. 



1 
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n 
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Fig. 2. Electron density versus chemical potential for Jh = 4, 6, and 10 (right to 
left), and J' = 0.02. MC results at j3 = 50, L — 20 for the spinless fermion model 
^ESF (circles) are compared with those for the DE model H (crosses). Error bars 
of the MC data are smaller than the symbols. The lines correspond to groundstate 
UHA. From Ref. [15]. 

Figure 2 shows the electron density versus the chemical potential fi for 
the Kondo model, the effective spinless fermions and groundstate UHA for 
Jh = 4, 6, 10. All three models give almost identical results. For very small 
H, the band is empty and J' leads to antiferromagnetism. At a critical \i\ 
depending on J c q = J', the filling jumps to n w 0.2 and the correlations 
become ferromagnetic. At a second critical /i!;, depending on J„g- = J' + 
1/(2 Jh), it becomes antiferromagnetic again and the density jumps from n w 
0.7 tonsil. 

These discontinuities mean that intermediate particle numbers are not 
stable in the grand canonical ensemble. They have been interpreted as phase 
separation. However, we will show in Sect. 3.3 that their cause lies rather in 
small ferromagnetic polarons. 

3.2 Finite Temperature UHA and Curie Temperature 

Although the uniform hopping approach replaces the fluctuating core spins 
by an average u, it reproduces not only the expectation value of the energy, 
but also its width with astonishing accuracy, even for higher temperatures. 
The results of UHA also remain valid upon the inclusion of n.n Coulomb 
repulsion [17]. All the while, the numerical effort is reduced from sampling 
over hundreds of thousands of core-spin configurations to scanning the sin- 
gle parameter u within the unit interval. The fact that the results remain 
valid with inclusion of the Coulomb repulsion indicates that UHA is a reli- 
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able approximative method that can safely be extended to more complicated 
situations. 
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Fig. 3. Curie temperature (dashed line) of the one-orbital DE model for a 16 3 cluster 
and t — 0.2 eV calculated in UHA. Circles and phases PM (paramagnetic metal), PI 
(paramagnetic insulator), FM (ferromagnetic metal), Fl (ferromagnetic insulator), 
and SCI (spin canted insulator) are experimental results for Lai-^Sr^MnOa [19], 
UHA results [17]. 



As observables can be evaluated for all temperatures once the density of 
states r(u) is known, it can be used to determine the Curie temperature for 
the three dimensional one-orbital model with Jh = oo and J' = 0. If one sets 
the only free parameter, namely the hopping strength t, which was also used 
as unit of energy, to t = 0.2eV, in accordance with experiments, one obtains 
the Curie temperature in reasonable agreement with experiment, see Fig. 3. 
In order to obtain the different low temperature phases besides FM and PM 
observed in experiments for low carrier concentrations, finite Jh and J' would 
be needed as well as two orbitals with Coulomb repulsion. 

3.3 Phase Separation versus Ferromagnetic Polarons 

The discontinuity of the filling as a function of the chemical potential, see Fig. 
2, is usually interpreted as phase separation [18], i.e. the system is expected 
to split into antifcrromagnctic domains with low and ferromagnetic domains 
with higher carrier concentration. Taking Coulomb interactions into account, 
PS has been argued to lead to either small [20] or large [21] (nano-scale) 
clusters, which have been the basis for a possible though controversial [22] 
explanation of CMR [2,23]. More thorough evaluation of the MC data for 
the transition near the filled lower Kondo band reveals however, that single 
hole ferromagnetic polarons are stabilized instead, even without any Coulomb 
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repulsion. All results presented in this section are for the effective spinless 
fermion model (7). 
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Fig. 4. Mean particle numbers Ne in a grand canonical MC simulation in 2 dimen- 
sions (L = 14 x 12, J H = 6, J' = 0.02, /3 = 50) as a function of MC time. One time 
step corresponds to 200 sweeps of the lattice, (a) pi = 1.26 > pi*: almost filled, (b) 
H = 1.19 ~ n*: polaron regime, (c) fi — 1.12 < pi*: FM regime. For visibility, only 
the first 350 time steps are shown. 



A first indication for ferromagnetic polarons is the behavior of the elec- 
tron density {n)s with the MC time near the critical chemical potential as 
depicted in Fig. 4, where (n)s is the thermodynamical expectation value for 
the filling given the core-spin configuration S sampled by the MC run. In 
the FM regime [Fig. 4(c)], the density fluctuates slightly with the fluctuat- 
ing core-spin configurations and takes non-integer values in accordance with 
standard results for free electrons. For the almost filled band [Fig. 4(a)], the 
density is (n)s = 1 (0 holes) for most spin configurations and occasionally, 
configurations occur, which contain exactly one, exactly two or exactly three 
holes. In between [Fig. 4(b)], the particle number fluctuates strongly, but as 
for (ri)s ~ 1, almost only integer fillings occur. While these integer fillings 
can hardly be understood in a PS scenario which is supposed to be a mix- 
ture of the low- and the high-density phases, it can easily be explained by 
independent polarons containing one hole each. As they are independent from 
each other, all have the same energy exactly balanced by the critical chemical 
potential pi* and their number therefore fluctuates strongly. 

To measure the size of the ferromagnetic domains, we use the dressed 
core-spin correlation function 

L-l 

S h {l) = — Y,n 1 l Si-S i+l (18) 
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that measures the correlations around a hole. The hole- density operator is 
related to the electron density via = 1 — rij. Equation (18) holds for open 
boundary conditions (employed in ID), the formula for periodic boundary 
conditions used in 2D is analogous. 




Fig. 5. Dressed spin-spin correlation function (18) from unbiased MC for j3 = 50, 
J' = 0.02, and Jh = 6. Left panel: 1 dimension: L = 50-site chain containing one 
(x), two (o), three (o), four (<), and five (>) holes. The dashed line is calculated 
within the simple polaron picture, while the solid line represents the generalized 
UHA result for a single polaron, see [24]. Right panel: 2 dimensions: 12 x 14 lat- 
tice with 1 (o), 6 (x) and 20 (A) holes. Continuous lines are data for the simple 
Polaron model (see Fig. 6, right panel): 1 Polaron (dotted), 6 (solid) and 20 Po- 
larons (dashed). The insets show the conventional spin-spin correlation function 

s (0 = t~i Y,i=l s * ■ S *+i- Left from Ref - I 24 !' ri § ht from Ref - I 25 ]- 

Figure 5 shows this dressed core-spin correlation. The ferromagnetic re- 
gions around the holes are small and their size does not grow with doping 
neither for the one- nor for the the two-dimensional case. This indicates that 
the introduction of more holes leads to more small FM polarons rather than 
to a growth of the existing ones. 

These observations lead to the development of a simple toy model for the 
FM polarons. Each polaron consist of a small (3-4 sites in ID, 5 sites in 2D) 
FM well, in which the hole can delocalize. These wells are embedded into an 
AFM background. For a schematic representation, see Fig. 6. The value of 
the critical chemical potential can be obtained from the toy model by simple 
energy considerations. It is simply the difference between the energy gained 
by the delocalized hole and the energy payed for the breaking of AFM bonds. 
In this simplest model, the hopping strength is given by uj = 1 for the FM 
regions and u a = for AFM bonds. The impact of thermal fluctuations of 
the core spins can be modeled by a generalization of UHA to inhomogeneous 
structures, where r(v,f,u a ) has to be determined. 

This was done for the one-dimensional model, for results and details of 
the algorithm see [24]. In order to compare the toy model to the MC data in 
2 dimensions, random deviations were added to the core spins, see [25]. The 
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Fig. 6. Toy model for FM polarons in one (left) and two (right) dimensions. Height 
represents hole density. In 1 D, a FM domain of Lf — 4 lattice sites is embedded in 
an AFM background, in 2D one spin is flipped from the perfect AFM. A single hole 
is localized in the FM domain giving rise to the depicted hole density (different from 
the schematic shape in Fig. 4 of Ref. [26]). Left from Ref. [24], right from Ref. [25]. 

principal effect of the fluctuations is a finite bandwidth for the AFM at half 
filling, their amount was therefore fitted to yield the same width for the AFM 
band as the MC data, see Fig. 8. 

Because the FM wells in which the hole can move are so small, they give 
rise to only a few well separated signals in the spectral density. Figure 7 shows 
the spectral density and the density of states for one, two and three holes in 
one dimension. On sees a broad band in the center which comes from holes 
moving in the imperfect AFM background. Separated from this central band 
by a (mirror) pseudogap are dispersionless states from the FM polarons at 
lo ~ ±1.5. The weights of these signals increase upon the introduction of 
more holes. 




LO U) LO 



Fig. 7. Spectral density for Nh = 1 hole, Nh = 2 holes Nh = 3 holes, corresponding 
to 1, 2 and 3 polarons in one dimension. Dashed lines: MC data, solid: generalized 
UHA. Parameters as in Fig. 4. From Ref. [24]. 

Figure 8 shows the spectral density for the two dimensional model with 6 
and 20 holes and compares it to the data for the toy model with added random 
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Fig. 8. Spectral density for J' = 0.02, (3 = 50, J H = 6 on a 12 x 14 lattice: left: 6 
holes (x w 0.035); right: 20 holes (a; w 0.12). Solid lines are the unbiased MC data, 
dashed lines the simplified polaron model. 

fluctuations. There is again a broad central band from the AFM featuring a 
mirror band due to the doubling of the unit cell and again the polaronic states 
separated by a pseudogap. There is strong correspondence of the unbiased MC 
data to the simplified model for both the one- and the two dimensional system. 

The pseudogap, which is also observed in experiments [27-30], can easily 
be explained by the few well separated eigenenergies of the holes trapped in 
the small polarons. In a phase separation scenario with larger FM regions, 
additional states would fill the gap between the AFM band and the polaron 
states, in contrast to experimental results and the MC data. 

3.4 Phase diagram in 2D 

Although < J' < 0.1 is by far the smallest parameter of the Hamiltonian, 
it has a considerable effect; especially at low carrier density, when the kinetic 
energy is small. While we observed nearly independent polarons for j3 — 50 and 
J' = 0.02, they tend to form larger clusters and eventually phase separation 
for decreasing J' and avoid each other for larger J'. The reason for this effect 
is the stabilization of the AFM background by J'. 

In order to determine the phase boundary between the polaronic regime 
and phase separation, we chose the filling, at which the nearest AFM signal 
(at the distance r 2 = 5) in the dressed core-spin correlation (18) became 
ferromagnetic. This criterion is somewhat arbitrary and the transition is not 
sharp, polarons rather coexisting with larger clusters. 

For large doping, J' — 0.05 destroys ferromagnetism and leads to the so 
called flux phase around half filling [11, 12,31,32]. In the phase diagram, one 
sees a small window in parameter space for phase separation, but for realistic 
parameters J' > 0.01, polarons dominate, especially at small doping. 
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Fig. 9. Electron number 7V e as a function of \x and J', and phase diagram for 
-0.5 < a* < 1-6, < J' < 0.05, < x < 0.6 (i. e. filled to 40% filled lower Kondo 
band),Jjj = 6, fi = 50 on a 14 x 12-lattice. "Pol.": polaronic regime, "Pol. /PS": 
mixture of both polarons and larger ferromagnetic clusters, "FM": ferromagnct, 
"AFM": antiferromagnet, "PM": regime without magnetic structure, "Flux": Flux 
phase. 

4 Summary 

An effective spinless fermion model (ESF) was derived from the FM Kondo 
lattice model (1) with classical core spins and for Jh S> t, leading to the ESF 
Hamiltonian (7). It allows the treatment of finite Jh with the same numerical 
effort as the Jh — > oo approximation, but gives better correspondence to the 
original model. A further simplification was achieved by treating the fluctuat- 
ing core spins by a uniform hopping strength (UHA) . With this much simpler 
model (13), we obtained the Curie temperature for the 3D model with one 
itinerant orbital in accordance with experimental values. 

By unbiased MC simulations for the ESF model (7) with non-degenerate 
conduction band in one and two dimensions, we found that ferromagnetic 
polarons are the reason for features previously attributed to phase separation. 
This polaronic behavior is enhanced by larger J' > 0.02. A phase diagram was 
obtained for the 2 dimensional case. 
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